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ABSTRACT 

We present an indirect dynamical evidence, based on the measure of the core to half mass 
radius ratio (r c /rh), that a significant fraction of globular clusters in our galaxy with an age 
greater than 10 half-mass relaxation times (t r h) host an intermediate mass black hole (IMBH). 
In fact, after a few t r h much of the memory about the details of initial conditions is erased and 
r c /rh is determined solely by the balance of the stellar encounters energy production in the 
core with the dissipation due to global expansion of the cluster. Here we compare the observed 
values of r c /r^ from a sample of 57 galactic globular cluster, selected to be dynamically old 
and not strongly influenced by tidal forces, with the theoretical expectation for the r c /r/, 
ratio based on analytical models and detailed numerical simulations. The simulations of the 
evolution of star clusters considered include combinations of single stars, primordial binaries 
and IMBH. For at least half of the clusters in our sample the observed r c /r^ ratio appears to 
be too large to be explained without invoking the presence of an IMBH at the center of the 
system. 



Key words: black hole physics - 
n-body simulations 



stellar dynamics — globular clusters: general — methods: 



1 INTRODUCTION 

The existence and detection for both stellar mass black holes 
as well as supermassive black holes are commonl y accepted 
and recognized in the ast ronomical community fad 1 19991 : 
iKormendv & Gebhardtl 1200 ll) . Intermediate mass black holes 
(IMBH), i.e. black holes with mass in the 10 3 - 10 4 Mq 
range, possibly formed by runaway mergers of massive stars 
jPortegies Zwart et alj|2004h are expected to fill in the mass range 
from stellar to supermassive BH. Indeed, evidence for the presence 
of IMBH in star cluster is accumulating, but so far there are only 
two globular clusters - M15 and Gl - that present observational 
evidence for the existence of IMBH jGebhardt et alJl2002L 120051 : 
Ivan der Marel et al.ll2003: iGerssen et alj|2003l) . The cases for the 
presence of IMBH are based on the challenging determination of 
the precise kinematics of the systems, that has to be compared with 
detailed dynamical models in order to infer the presence of the cen- 
tral BH. At the level of velocity dispersion and full line of sight 
velocity dispersion profiles, the direct dynamical influence of an 
IMBH is limited within the sphere of influence of the BH, hard 
to resolve observation ally (but see the recent hig h-resolution HST 
survey carried out by iNovola & Gebhardl 2006). In addition, the 
intrinsic three dimensional density (and velocity dispersion) cusp 
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around the central BH is smeared out when the data are projected 
along the line of sight. 

However, even if the direct dynamical signature of a central 
BH is hard to detect observationally, the presence of an IMBH 
at the center of a collisional stellar systems has profound conse- 
quences on the global dynamics. In particular the visual appear- 
ance of a globular cl usters hosting an IMBH is that of a system 
with a sizeable core ( B aumgardt et all 120051 : iTrenti et al.1 [2006b ) 
and not that of very concentrated cluster, as could be naively in- 
ferred from the existen ce of the Bahcall-Wolf 1/r 1 ' 7 density cusp 
(Bahc all&Woiill 19761) within the sphere of influence of the BH. 
This is because the presence of a central BH and of its compact den- 
sity cusps generates, through stellar encounters, a sufficient amount 
of energy to halt the core collapse, if the initial conditions start 
with a wide core, or to fuel a core expansion if t he initi al config- 
uration is centrally concentrated. In ITrenti etHll d2006bh we have 
shown that star clusters with a central IMBH starting from a va- 
riety of initial stellar density profiles (i.e. King models with con- 
centration parameter Wo = 3, 5, 7, 11 as well as Plummer models) 
all evolve toward a common value of the core to half mass radius 
ratio (r c /rh ~ 0.3) on a timescale of a few initial half-mass re- 
laxation times. This is indeed not surprisingly: a similar behavior 
is observed when primordial binaries are the responsible driver for 
energy production in the center of the cluster dTrenti. Heggie & Hull 
2006), in agreement with the t heoretical expectation formulated 
bv lVesperini&ChernofJ < fl994n . As binaries alone are less effi- 
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cient at producing energy than an IMBH, in this case we have 
Tc/th < 0.06 (where the upper limit is reached in presence of a 
number density of at least 50% binaries in the core) for a typical 
globular cluster. Only when no or little binaries are present there 
can be a deep core collapse , ending when r c /rh < 0.02 (see, e.g. 
iHeggie. Trenti & Hutll2006h . 

Given this theoretical framework, the observed core to half 
mass radius ratio appears indeed a powerful indicator of the dy- 
namical configuration of the core of an "old" globular cluster, i.e. 
a globular cluster that has evolved for about 5 to 10 half-mass re- 
laxation times. r c /rh is in addition a relatively easy photometric 
measurement, that can be applied even to extra-galactic globular 
clusters, without worrying about the need of a very accurate mea- 
sure. In fact, the combination of numerical simulations and theoret- 
ical modeling predicts that this ratio increases by about a factor 3 
moving from clusters with single stars only to clusters with a size- 
able population of binaries; when an IMBH is present r c /rh is 4 to 
5 times larger than in the case when only binaries are present and 
12 to 15 times larger than when the cluster contains single stars 
only. 

In this letter we present a preliminary application of the use of 
r c /rh as an indirect indicator for the presence of an IMBH in the 
cores of 57 galactic globular clusters that satisfy a conservative dy- 
namical age criterion, i.e. their half-mass relaxation time is shorter 
than 10 9 j/r. In particular, we show that the observed core to half 
mass radius ratio in more than half of the clusters in the sample is 
so large that the only general dynamical explanation can be given in 
terms of the presence of a central BH. The paper is organized as fol- 
lows: in Sec.|2]we summarize the current numerical and theoretical 
understanding of the evolution of r c /rh, in Sec.[3]we present the 
sample of galactic globular clusters that we use in this analysis and 
we discuss their properties in terms of our theoretical expectations. 
Finally we conclude in Sec. [4] 



2 GLOBAL EVOLUTION OF A CLUSTER WITH 
BINARIES AND IMBH 

As discussed in quantitative detail in a series of paper devoted to the 
dynamic al evolution of star clusters with primordial binaries and 
IMBHs dHeggie. Trenti & Hu j l200d : ITrenti. Heggie & Huj|2006l : 
iTrenti et aUl2006bh . the long term evolution of a globular cluster 
is driven by the heat flow generated in the core of t he system to- 
ward the halo (see also IVesperini & Chernofdll994l) . Three main 
processes can be identified at the base of this heat flow, that reg- 
ulates the evolution of r c /rh (shown in Fig. Q] for representative 
simulations of our numerical simulations program). (1) When only 
single stars are present the core behaves as a self-gravitatin g sys- 
tems with negative specific heat (e.g. see lHeggie & Hutll2003l) and 
a thermal collapse happens leading to a core contraction. The core 
contraction proceeds for about 10 — 15 half mass relaxation times 
for equal mass particles, while it is significantly faste r when a mass 
spectrum is considered (Chernoff & Weinberg 1990). The collapse 
lasts until the central density is so high that a few binaries are dy- 
namically formed by three body encounters. Given the high cen- 
tral density these binaries interact efficiently with single stars and 
become progressively tighter and tighter providing enough energy 
to halt the collapse and create a core bounce. When the number 
of stars in the system is large enough (N > 10 4 ) binary activity 
may cause a temporary temperature inversion in the core and subse- 
quent core expansion. The proc ess is eventually repeated in a series 
of "gravothermal oscillations" dSugimoto & Bettwieserll983t) with 



the core radius oscillating around a value roughly 2% of the half 
mass radius. (2) When a sufficient fraction (i.e. greater than a few 
percent) of stars in the cluster has initially a companion with typical 
separations below 10 AU (that is "hard binaries" that have a binding 
energy of at least a few times the total energy of the cluster), energy 
generation due to existing binaries is much more efficient than that 
due to dynamically formed pairs. Here the equilibrium size of the 
core is larger than when only single stars are present. The system 
may even show a core radius expansion if the initial conditions are 
too concentrated. The precise value of the equilibrium core radius 
has a moderate dependence on the number of stars in the systems 
and on the number density of binaries in the core: r c ~ 5% for a 
typical glo bular cluster (N = 3 ■ 10 5 ) with about 40% of binaries in 
the core dHeggie Trenti & Huj|2006l: ITrenti. Heggie & Hutll2006l: 
iFregeau & Rasiol |2006t IVesperini & Chernofj 1 19941) . This quasi- 
equilibrium evolutionary phase lasts until there is a sufficient num- 
ber of binaries in the core, that is until the binary fraction in the 
cluster is reduced at a few percent level. In fact, binaries in the outer 
part of the cluster tend to sink toward the core for mass segregation, 
so the core binary fraction is replenished until the total binary frac- 
tion is almost exhausted. (3) When an IMBH with mass of the order 
of 1% of the total mass of the cluster is present, a yet more efficient 
energy production channel is available in the center of the system 
via dynamical inter actions within the influence sphe re of the BH 
(Trenti et al. 2006B see also lBaumgardt et alj|2004al lbl). The equi- 
librium value of the core radius in this case is significantly larger: 
r c /rh ~ 0.3 independently from the initial binary fraction. Even 
in this case there is an expansion of the core radius on a relaxation 
timescale if the initial conditions are too concentrated. 

This general picture has been obtained with somewhat ideal- 
ized N-body simulations that employed single particles only and 
did not take into account stellar evolution. Therefore it is impor- 
tant to assess if the results that we obtain can be biased by our 
assumptions. Modest deviations from our numerical expectations 
will not impair the usefulness of r c /rh as diagnostic tool to infer 
the presence of IMBH, as r c /r^ is expected in this case to be larger 
by a factor 4 to 5 with respect to the case where only binaries are 
present. 

To ensure that no major bias is present we examine critically 
the idealizations in the numerical runs: 

(i) Mass spectrum. The choice to use equal mass stars is prob- 
ably the most significant simplification introduced in our simula- 
tions. When a realistic mass spectrum is present, massive stars seg- 
regate on a relaxation timescale toward the core of the system and 
speed up the core collapse by ab out one order of magnitude, i.e. 
the collapse takes about 1 — 2 t r h (Chernoff & Weinberg 1990). In 
principle, the equilibrium value of the core radius could be changed 
significantly. However simulations of the dynamics of star clus- 
ters with primordial binaries that include a realistic mass spec- 
trum have been recentl y performed with a Monte Carlo code by 
IFregeau & Rasiol d2006h and their results on r c /rh are fully con- 
sistent with our more idealized runs. Similarly, realistic simula- 
tions of globular clusters with an IMBH that include a mass spec- 
trum a nd stellar evolution (but w ithout binaries) have been carried 
out by Baumgar dt et al.l d2004bl) and again the presence of a large 
core radius is confirmed. The anal ysis of the final configura tion at 
t — YlGyr of two largest runs bv lBaumgardt et al.l d2004bh - with 
N = 131072 - leads to r c /r h « 0.27 for the run with a lKroupal 
d200ll) IMF truncated at 100 M & and r c /r h » 0.35 when the IMF 
is truncated at 30 Mq. 

(ii) Stellar evolution. Our runs do not take into account the ef- 
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fects of stellar evolution, that limits the lifetime of massive stars 
and may lead to the coalescence of tight binaries. These effects 
may possib ly lead to a more r apid depletion of the binary popu- 
lation (see llvanova et alj|2005h . so that a star cluster would burn 
out its binary reservoir earlier than expected from our simulations. 
If this is the case, then the core undergoes a deep core contrac- 
tion, like in the case where only single stars are present. The di- 
rection of the evolution is fortunately in the right direction to avoid 
an observational bias when we are interested, as in this paper, to 
focus on large values of the core t o half mass radius. For s imu- 
lations with an IMBH, the study of iBaumgardt et alj d2004bl) (see 
point above), that included stellar evolution, guarantees instead that 
the bias due to stellar e volution is not important (see also Fig. 1 in 
IBaumgardt et alj2004bl . where the Lagrangian radii, that is the radii 
enclosing a fixed fraction of the total mass of the system, expand 
steadily starting from a Wo = 7 model). 

(iii) Number of particles. Our direct simulations have been lim- 
ited by CPU speed to employ only up to ~ 20000 particles. 
The typical number of stars in globular clusters is larger by at 
least one order of magnitude. Given our choice to use equal mass 
particles and no stellar evolution, our results are scale-free and 
can be adapted to any physical value for the mass and radius 
scales . From theoretical considerations by Vesperini & Chernoff 
d 19941) when primordial binaries are present r c /rh is expected to 
have a modest dependence on the number of particles N, that 
is r c /rh oc l/log(0.11N). Simulations with N from 512 to 
16384 i ndeed verify that the scaling is app roximately consistent 
with the IVesperini & Chernofdll994l model (He ggie. Trenti & Hutl 
l2006h . If anything r c /r h (N) seems to decrease sligh tly /aster than 
expected. In addition the simulations by iFregeau & R asio 
employ a realistic number of particles and t heir core radii are con- 
sistent with our results extrapolated using the lVesperini & Chernofj 
d 19941) model. When a IMBH is present we do not expect a sig- 
nificant scaling of r c /rh with N as confirmed by our a nalysis of 
the two N = 131072 runs of IBaumgardt et alj f2004bT) . Even if 
a log ./V scaling were to be assumed despi te the indication from 
the simulations bv lBaumgardt et all d2004bl) . the core radius value 
extrapolated from our runs with a IMBH would still be such that 
fc/Vh > 0.15 for a realistic number of particles. 

(iv) Tidal field. The runs with an IMBH in lTrenti et all j2006b) 
do not take into account the presence of a tidal field. This is how- 
ever not expected to introduce a significant bias: r c /rh in runs with 
primordial binaries varies at t he 10% level with or witho ut inclu- 
sion of the galactic tidal field dTrenti. Heggie & Hutll2006h . 

(v) Spherical symmetry. Our initial conditions are all spherically 
symmetric. We do not expect this to be a problem as globular clus- 
ters are typically very close to spherical symmetry. In addition, the 
precise details of the starting configuration are not important on a 
collisional timescale. To be sure, in the analysis of the next sec- 
tion we analyze a subsample of the clusters selected by excluding 
objects that departs from spherical symmetry, obtaining the same 
rc/rh distribution as in the parent sample. 

(vi) Rotation. The simulations that we consider have zero net an- 
gular momentum, i.e. no rotating model is studied. Observational 
evidence of rotation in globulars is weak and t he velocity disper- 
sion t ensor displays a mild anisotropy at most dMevlan & Heggie] 
1 1997b . Therefore it would be very surprinsing if the comparison 
with observations turns out to be strongly biased by not including 
rotation in simulations. 



2.1 Core radius definition 

Finally, before moving on to compare our theoretical expectations 
for r c /rh with the data, we have to consider a yet different source 
of uncertainty. Numerical simulations give us mass-defined core 
and half-mass radii, while luminosity-based quantities are derived 
from observations. A discrepancy can therefore arise (1) if the lu- 
minosity profile does not track accurately the mass profile and/or 
(2) if the core radius is differently defined in the two cases. 

For the first issue we have to consider that the luminosity of a 
globular cluster is dominated by stars in the giant branch, that are 
more massive than the average star of an old globular cluster. The 
most luminous stars are more centrally concentrated due to mass 
segregation than the average, so if a bias is introduced, it goes in 
the direction of reducing the observe d core radius with respect to 
a mass weighted definition. However lChernoff & Weinberg! dl990l) 
derive only modest color gradients due to this effect. In any case it 
does not affect the interpretation of large observed cores. 

The second issue, that is a difference between the observa- 
tional and theoretical core radius definition is also not expected 
to be a major problem. In our simulations we adopted the den- 
sity w eighted definition of the core radius from Case rtano & Hutl 
dl985l) . Eq. IV.3: 

_ (W\p)m = E t npirrii 

C ~ (P)M Pi m i 

where the sum is carried over all the stars in the simulation, n 
is the distance of the i-th star from the center of the system, rrn 
its mass and pi the stellar density computed at the star position. 
This definition for r c is closely aligned to the standard observa- 
tional practice of defining the core radius as the radius r M where 
the luminosity surface density has dropped to half its central value 
dCasertano & Hutlll985t) . 



3 DATA SAMPLE 

As a preliminary application of the r c /rh diagnostic to infer the 
presence of an IMBH, we consider the c ompilation of galactic glob- 
ular clusters properties by lHarrisId 19961) , revised in Feb ruary 2003. 
For each globular cluster in the tables o f lHarrisI < fl99^> we extract 
the following quantities: (1) core radius (r c ); (2) half-mass radius 
(rh); (3) tidal radius (r t )\ (4) half-mass relaxation time (i r h); (5) 
ellipticity (e); (6) cluster luminosity - i.e. the absolute visual mag- 
nitude - (M vt ). 

We then proceed to build a uniform and homogeneous data 
set that includes only old globular clusters to ensure that these sys- 
tems are well relaxed by two body encounters. As a first step we 
therefore exclude from the analysis: 

(i) All the globular clusters for which one of these values is not 
quoted (except for e); 

(ii) All the globular clusters with t r h > W 9 yr. Assuming a typ- 
ical globular cluster age of 10 10 yr, this ensures that all the objects 
in the sample are at least 10 half mass relaxation times old. Most of 
the clusters excluded from the sample fail this test; 

(iii) All the cluster that are under the influence of a strong tidal 
field effects, i.e. those for which r t /rh < 4. This last selection 
criterion is mainly a precaution to exclude systems where r c /rh 
may s tart fluctuating as rt approaches m (see lTren ti. Heggie ^; Hutl 
2006). Only four clusters that pass the previous selections are ex- 
cluded due to a small tidal radius. 

After applying our selection criteria we are left with a sample 
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Figure 1. Evolution of r c /rh i n a series of runs with e qual m a ss stars from 
Heggie , Trenti & Hull 120061) ; iTrenti, Heggie & Hud 120061) ; iTrenti et alj 
12006b) starting from a variety of initial conditions that include single stars 
only (red line: Plummer model, N=8192), primordial binaries with galactic 
tidal field (green line: King model Wo = 7 & 10% binaries, Wo = 11 
& 20% binaries, N = 16384) and with IMBH (black line: m hh = 
lA%M to t & 10% binaries, Plummer model, N=8192; m bh = 3%M tot , 
no binaries, King model Wo = 11, N=4096). When a realistic mass spec- 
trum is considere d for single stars runs the core collapse is much faster 
(rj 2 t rh , see lChernoff & Weinberg! 19901) . 

of 57 old globular clusters, whose r c /rh distribution is plotted in 
Fig-El The plot is extremely interesting: it is immediately apparent 
that only a modest fraction of the globular clusters is character- 
ized by r c /rh < 0.1, that would be consistent with the core size 
expected when the system is populated only by single and binary 
stars. The majority of the objects in the sample has a large core 
(fc/rh > 0.2) and r c /rh > 0.5 is not infrequent. The average 
value of the ratio is (r c /rh) = 0.31. The median of the distri- 
bution is at r c /rh = 0.28. From the dynamical considerations of 
Sec. [2] only the presence of an IMBH appears to be consistent with 
such large core radius values. 

In Sec. [2] we have seen that, when the core to half mass ra- 
dius ratio is set by the burning of binaries there is a modest, log- 
arithmic N- dependence on this ratio, th at is reasonably well cap- 
tured by the lVesperini & Chernoff ( 1994) m odel, especially in p res- 
ence of a galactic tidal field (see ITrenti. Heggie & Hutl [20061) . To 
take into account this factor in the analysis we assume (1) that the 
total luminosity is a proxy for the total mass of the system and 
(2) that the average stellar mass in each cluster is constant and 
equal to 0.7Mq . The mass-lumino sity relaxation has been ch ecked 
and calibrated using the data from lGnedin & OstrikeJ Jl997h . We 
used log 10 (M to t/M ) = A + b ■ M V t with A = 2.5 and 
b — —0.38236. We have then corrected the ratio r c jru by applying 
a correction factor £ 

log (0.11/0.7 M tot /M Q ) 
? ~ log (0.11 AW) 
so as to obtain a ratio f c /r~h — (,r c /rh equivalent to that of a clus- 
ter with N re f = 3 • 10 , a standard value in numerical simulations 



of globular clusters. The results shown in Fig. [2] do not change 
significantly when this correction is applied: the average value is 
{f c /r~h) — 0.27 and the median value of the distribution is 0.22, 
only marginally smaller than without the correction. The presence 
of an IMBH is therefore still needed to explain the large value of 
the core radius for at least half the objects in the sample. 

Could this conclusion be biased by the fact that the observed 
globular clusters are for some (improbable) reasons not yet well 
relaxed, so that the equilibrium value for r c /rh has not yet been 
reached after ~ 10 t T ht To address this point we enhance our se- 
lection algorithm to select only objects that are at least 20 t r h old, 
i.e. with t r h < 5 • 10 8 yr. This reduces our sample to 25 clus- 
ters whose average core radius is still (r c /rh) = 0.36 (marginally 
larger than the whole sample value), while the median is at r c /rh = 
0.21 (marginally sma l ler th an the whole sample analysis). If the 
Vesperini & Chernof3 l l 19941) Zo<?(0.11./V) correction is applied we 
have {f c /r~h) = 0.28 and median value of f c /r~h = 0.17. The net 
effect of the correction is to reduce the mean and the median val- 
ues, as by requiring t r h ^ 5 ■ 10 s yr we are biased toward selecting 
smaller clusters. Even for this reduced sample more than half of the 
objects have f c /r~h > 0.15, so the need for a central IMBH remains 
strong. 

As an additional test to ensure absence of observational biases, 
we have further refined the reduced sample with t r h < 5 • 10 8 yr to 
exclude globular clusters that depart from spherical symmetry, i.e. 
those clusters with a quoted e > 0.1, where the ellipticity ratio e 
is defined in terms of the axis ratio (e = 1 — a/6). This leaves 20 
objects in the sample and does not change the results: (r c /rh) = 
0.38 (median 0.23), while (f c /r h ) = 0.29 (median 0.20). 

All the progressively selective cuts that we applied to the 
galactic globular cluster system provide a consistent picture where 
about half of the objects have a significantly large core radius, i.e. 
Tc/rh > 0.2. If an IMBH turns out to be absent in these sys- 
tems, then a different dynamical explanation for their large core 
radii must be found. If this is indeed the case for the whole sam- 
ple then a major rethinking of our knowledge of the dynamics of 
globular clusters is needed. Nevertheless a few individual objects 
in the sample may turn out to have a large core radius for pecu- 
liar reasons, especially if they have been strongly perturbed and are 
therefore not as dynamically old as their relaxation time would im- 
ply. One example of a peculiar object that passes all our selection 
criteria is the cluster with the largest r c /rh ratio, Pal 13, which is 
considered a very unusual object, either due a strong tidal heating 
during the last perigalacti con passage or d ue to the presence of a 
sizeable dark matter halo dCote et alj2 002). 



4 CONCLUSIONS 

In this paper we propose the use of the observed ratio of the 
core to half mass radius as a powerful indirect dynamical in- 
dicator for the presence of an IMBH at the center of old stel- 
lar cluster. A numb er of theoretical and numerical investigations 
combi ned together |Ve sperini & Chernoff 1994; Baumgar dt et al.l 
2004allbl: iFregeau & Rasiol 12006): iHeggie. Trenti & Hull 120061: 
Trenti, Heggie & Hut 2006; ITrenti et ai1l2006bf) strongly support 
the idea that after about 5 to 10 half-mass relaxation times the 
value of r c /rh in a globular cluster assumes significantly differ- 
ent values depending only on whether the core contains single stars 
{fc/rh ~ 0.02), a large fraction of hard binaries (r c /rh ~ 0.05) 
or an IMBH (r c /rh ~ 0.3). The details of the initial conditions, 
such as the initial density profile, do not affect this quantity. By an- 
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Figure 2. Main panel: distribution of the r c /r^ ratio in the main sample 
of 57 old galactic globular clusters (selected to have t r ^ < 10 9 yr). Small 
panel: £ ■ r c /r h (corrected to rescale it to a reference number of particles 
N re f = 3 ■ 10 5 ) for a sub-sample (20 objects) of the main sample, selected 
by imposing t r f l < 5 ■ 10 8 yr and e < 0.1. Both panels show a significant 
number of clusters with Tc/r^ > 0.2. 



alyzing the distribution of r c /rh in a sample of 57 galactic globular 
clusters, selected to ensure that they are well collisionally relaxed, 
we conclude that there is a strong indirect evidence for the presence 
of an IMBH in at least half of the objects in the sample. 

A globular cluster that hosts an IMBH represents a laboratory 
where dynamical interactions between hierarchical systems take 
place frequently, leading in particular to the formation of triples 
(and even quadruples) systems and to the ejections of a number of 
high velocity stars, with som e able to reach ejec tion velocities up 
to several hundreds km/s (see lTrenti et alj2006bl) . 

Intriguingly, some objects in the sample have a core radius 
that could possibly be too large (r c /rh > 0.5) to be explained by 
the presence of a single IMBH. It could therefore be possible that 
some of these globular clusters host a binary IM BH. In this case, 
howe ver, the heating may be so efficient (e.g. see lYu & Tremainel 
120031 for an estimate of the interaction rate of a binary BH with 
single stars) that it is not clear whether the system is able to sur- 
vive for even a few relaxation times. Therefore an accurate numeri- 
cal modeling is critically required before any conclusion about this 
speculation can be drawn. 

Our conclusions on the presence of singles IMBHs appear in- 
stead to be robust. We have discussed a number of possible biases 
in the comparison of numerical simulations with the observations, 
but while some of the idealizations introduced in the modeling may 
induce limited changes in r c /rh, we could not identify a possible 
major bias. An error of the order at least 300% would be required 
to be able to explain the observed distribution of r c /rh without in- 
voking the systematic presence of an IMBH in half of the objects 
in the sample. 

Clearly more detailed numerical simulations will be required 
to evaluate and confirm the presence of IMBH in specific clusters of 
the sample. In particular it would be extremely interesting the con- 



version of N-body snapshots into synthetic observations that could 
therefore be directly analyzed like images acquired by a telescope, 
avoiding an indirect comparison between theoretical and observed 
quantities. 
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